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Q ' Abstract 

Q^ . Non-transparent models of multipole expansion model and two point- 

mass model are analyzed from the catastrophe theory. Singularity be- 
On' haviours of 2 n -pole moments are discussed. We apply these models to triple 
quasar PG1115+080 and compare with the typical transparent model, soft- 
ened power law spheroids. Multipole expansion model gives the best fit 
52 | among them. 
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1 Introduction 

Gravitational lensing is and will be playing very important roles in astrophysics 
1 1 1 . It covers the large region in scale from the brown dwarf to the supercluster 
of galaxies. Microlensing uncovers the dark matter in the halo and the planetary 
system in the neighbourhood of the Sun || . In this letter we deal with a special 
macrolensing, that is, the multiple quasar supposed to be caused by a single 
galaxy in the foreground. It enables us to reveal the structure of the lensing galaxy 
without the ambiguity of surrounding galaxies. It also gives us the information of 
the large scale structure of the universe through the angular diameter distance. 

One of the most popular models of lensing galaxy is the spheroid with softened 
power law behaviours [|]]. The physical implication of this model is very clear 
and it is easy to compare with observations, such as luminosity distribution and 
ellipticity, etc. However these merits are also the demerits of this model since it 
is too restrictive to model shape. We are not sure whether the lensing galaxy is 
elliptical or spiral. It may be even complex system composed of several galaxies. 
This model is also analytically non-tractable. On the other hand multipole ex- 
pansion model (hereafter MPE model) is more general than the spheroid model. 
MPE model is applicable to any complex system as far as the lensing objects are 
compact with respect to impact parameters. Also it is analytically tractable. 

Using it, we are able to estimate the structure of the lensing object systemat- 
ically. That is, the more we incorporate multipole components, the more we can 
know the detailed structure of the lensing object. 

From the observational side, gravitational lensing effect seems to be very 
promissing. Deep space survey by HST is now working and Sloan Digitized Sky 
Survey will work soon ||. However the present status of observation is not still 
sufficient to uncover the affirmative results. For instance, parameter fitting is 
very sensitive to the image positions whose observation cannot be precise in 
many cases. We need more detailed and more precise data to be adjusted. In 
such a situation it is also very useful to study the same lensing phenomena by 
several models, which reveal the common and therefore true features of lensing 
objects. 

In the previous paper we discussed about the relationship between multi- 
pole moments and caustic singularities, and applied it to PG1115+080. In this 
treatise it is also shown that complex coordinates are very useful. 

The merit of MPE model is that we can probe into the detailed structure 
systematically by steping up higher order of multipoles. An arbitrarily complex 
but compact lensing object is analytically tractable as far as each 2 n -pole mo- 
ment is discussed separately. So in this article we study the structure of 2 n -pole 
term systematically. If we sum up all contributions, its treatment goes beyond 
analytical survey and is devoted to the numerical calculations. 

Two point-mass lens is situated at the antipodal point of multipole expansion 
model. 



Two point-mass model consists of a series of multipole moments. However 
it allows analytical survey as a whole though a straightforward application of 
catastrophe theory still needs numerical calculations. In order to understand the 
essential physical implication of any event, analytical survey is indispensable. 

We apply the above mentioned models to PG1115+080 and discuss their 
common and discrepant characters among them. 

This paper is organized as follows. In section two we review the metamor- 
phoses of caustic for later use whose detailed arguments are given in our previous 
work [|7| . This metamorphoses are applied to MPE model in section three and to 
the two point-mass model in section five. Section four is devoted to the analysis 
of the detailed structure of 2 n -pole contribution in MPE model. MPE model is 
applied to the multiple lens system of PG1115+080 in section six. The obtained 
numerical results are compared with other models, especially with elliptical lens 
model in section seven. Section eight is devoted to discussions. 

2 Metamorphoses of the Caustics 

For later studies, we summarize the metamorphoses of the caustics. This is 
mainly based on our previous work fTj. 

Using complex variables, singularities are classified in compact form. 

Critical line is 



D = -4det(0 a/3 ) 

= 4(|0 22 -| 2 -|0 2 /)=O , (a,/3 = z,z) (2.1) 

where <fi is Fermat potential and z is complex representation of image position, 
z = Xi + ix 2 - z is complex conjugate of z. Subscripts a, [3 on denotes the 
derivative d a dp(f) and so on. 

Here we use the Poisson equation 

<j>zz = "(I - «) • 

K is the normalized surface mass density and is equal to zero on the light path 
for non-transparent model which we will consider in this article. 
Then the critical line is reduced to 

<p zz ± ~ = (2.2) 

by the diagonalization in real coordinates. + (— ) corresponds to <f>u = 

(022 = 0). 

Cuspoid sequences are characterized by 



± Tdr = (2.5) 



L n <p z = (n = l,2,3) . (2.3) 

L is the derivative operator along the critical line, which is given by 

L = Ud z Dd z -d z Dd z ) . (2.4) 

o 

Their explicit forms are 

d 3 <j) 

dz 3 dz 3 
for the cusp (n = 1) singularity, 

d 4 (j) d 4 ^ d 3 (j) 2 

( ^ + ^ )±6 fe } =° (2 - 6) 

for the swallowtail (n = 2) and 

for the butterfly (n — 3). 

It should be noticed that the higher cuspoids must also satisfy the conditions 
of the lower ones. For instance, butterfly must satisfy Eqs. Q2.2|) , ( [2.5|) and (|2.6| ) 
as well as Eq. (|2T7|). 

Umbilic is 

<Pa(3 = . (2.8) 

As is easily checked this singularity is not realized for the non-transparent 
model. 

Beak-to-beak and lips are characterized by 

d z D = , (2.9) 

a) Lips A = -det£> Q/3 > (2.10) 

b) Beak-to-beak A < . (2.11) 

The explicit form of Eq. (|2.9|) is 

d 3 <p 



dz 3 
A takes the form 



. (2.12) 



d 4 6 
A = -4|^|<0 (2.13) 

for non-transparent model. Therefore Eq.( |2.12j ) leads to beak-to-beak. 



3 Multipole Expansion (MPE) Model 



Fermat's potential is 
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Using the expansion form of logarithmic function, 



(3-1) 
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we obtain the multipole expansion form of Fermat's potential, 
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Thus MPE model becomes a non-transparent model from the convergence 
requirement of Taylor expansion. This property seems to have interrupted pop- 
ularity in gravitational lensing. However this model is very useful and applicable 
to many lensing events and deserves further investigations as will be discussed. 

We can transform the potential into, 
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For this potential, critical condition and cusp condition become 
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(3.7) 

(3.8) 

(3.9) 

(3.10) 



where a n = (n + 2)/2 and /3„ = n(n + 1) . Im (Re) denotes the imaginary 
(real) part. 



4 The Properties of 2 n - Pole 

Unfortunately, if we consider the terms up to quadrupole moment, MPE model 
is already beyond analytical calculation. However, analytical survey is desirable 
for the essential physical interpretation. So we consider the case of 2 n -pole only, 
which allows analytical survey. 
Introducing the variables, 

— = 4Xe me , Y = VXcosn6, 

Eqs.( fT8D , (|3.9|) and (|3.10|) are rewritten into, 

\Q n \^X 1+ i=X + 2[3 n Y + f3 2 n , (4.1) 

Im(cusp) = smn9{(2a n - 3)X 2 + 2(3 n (a 2 n - 3)XY 

+f3 2 n (3a 2 n - 6a n + 1)X - (3 2 n (2Y + a n /3 n ) 2 ] = , (4.2) 

Re(cusp) = X 3 + (2a n + 3)(3 n X 2 Y + (3 2 {2(a 2 n + 3)XY 2 - (a 2 n - Qa n + 3)X 2 } 
+3[3 3 n (a 2 n + 2a n - l)XY + a n /#(3a„ - 2)X + P 3 n Y(2Y + a n f3 n ) 2 > , (4.3) 

X > Y 2 . (4.4) 

The last condition follows from the just definition of X and Y. As for the 
solution of (X, y), there are 2n solutions for X > Y 2 on lens plane and n solutions 
for X = Y 2 . 



At first, we examine the cusp conditions, Eqs. Q4.2|) and ( |4.3|) , and Eq. 
When sinn^ = 0, i.e., X = Y 2 , Eq.(^3]) becomes 

Re(cusp) = Y(Y + a n /3 n ) 2 (Y + (3 n ) 3 > . 

Therefore a part of region that satisfies the cusp conditions is X = Y 2 for Y > or 
Y < -/3 n (Y ^ -a n (3 n ). Further at the points, (X, Y) = ((3 2 , -j3 n ) , (c^, -a n /3 n ), 
the beak-to-beak condition is held. 
When sin n9 ^ 0, Eq.( [4.2| ) becomes 



(2a n -3)X 2 + 2[3 n (a 2 n -3)XY + [3 2 n (3a 2 n -6a n + l)X-f3 2 (2Y + a n [3 n ) 2 = . (4.5) 

After some calculations, we can obtain the region that satisfies Eq.(fO|) and 
Eq.i 



X > , -^- > Y > -A, . (4.6) 

2a n - 3 

We call the characteristic points, Pi, P2 and P3 defined by 



Pi = (13^ -p n ), P 2 = {a 2 Jl -a n (3 n 



P 



3 



n 



(2a n -3) 2 ' 2a n -3' 
The region in which the cusp conditions is satisfied is depicted in Fig.l for 



Fig.l 



For the other values of n, topological structure does not change. 
Next we examine the qualitative behaviour of the critical condition, Eq. fl4.lp . 
Firstly, irrespective of the value of \Q n \, the curve Eq. fl4.1|) passes through the 
definite point, 

(0 , - %) , 



and passes above Pi, 



tfl+4/n 
Y(X = (3 2 n ) = -A. + |Q n | 4 /n^n__ > _p n . 

For X > 0, the curve has a locally minimum point and Y tends to oo at X = oo. 
Furthermore, for fixed X, the value of Y increases as \Q n \ does. 

From the above arguments, the number of the solutions of the critical condi- 
tion and the cusp conditions changes when the curve Eq. (|4.1| ) 

case (0) is tangential to X = Y 2 for Y < —f3 n or, 

case (1) passes through P 3 , or 

case (2) is tangential to Im(cusp) = for Eq. fl4.6D - 

Each case is depicted in Fig. 2 for n = 2,4. 

Fig.2(a) 



Fig.2(b) 



For the case(O), the value of \Q n \ is 

i Q -i = i^' 5 B(B+ i)(n + 2) (;rT2)'- (4J) 

and the tangential point is P%. Therefore beak-to-beak singularity appears. 
For the case(l), the value of \Q n \ is 



\Qn\ = |Q?>| = ^i( 3 ^^) f = T'\n - l)|QJp| . (4.8) 



As known from graphical consideration, then cusps split into more cusps or 
merge to less cusps. Therefore butterfly singularity appears. 

For the case(2), only when n < 3, the lines of Eqs. ( fOP and (|4.2|) are tangential 
to each other in the region Eq. (|4.6| ) and the value of \Q n \ is 



\Qn\ = \Q^\ 



n 



2 + 9n + 16 + nVn 2 + Wn + 17 



2n(n + l)(n + 2) 
n(n 



, ..,.. + 1 + Vn 2 + 10n + 17) -,n/2 

x 1 4(^T^ J • (49) 

And then cusps merge to fold, i.e., swallowtail singularity appears. 
We surmmarize the result: 

(I) for n = 1, P 3 = (oo , oo) 



(i) |Q«I < !-»„ 

There are three solutions on J = 7 2 line and two solutions for X > 

y 2 , so 3 + 2x2 = 7 cusps on lens plane. 

(ii) |Q(°>|<|Qn|<|Q( 2) | 

After beak-to-beak singularity appears, there is a solution on X = Y 2 
line and two solutions for X > Y 2 , so 1 + 2 x 2 = 5 cusps on lens 
plane. 

(hi) |Qi 2 )| < \Q n \ 

After swallowtail singularity appears, there is a solution on X = Y 2 
line and no solution for X > Y 2 , so one cusp on lens plane. 

(II) for n = 2, 3 

« IQn|<|Q( 0) | 

There are three solutions on X = Y 2 line and one solution for X > Y 2 , 
so 3n + 2n = 5n cusps on lens plane. 

(ii) \Q®\<\Qn\<\Q$>\ 

After beak-to-beak singularity appears, there is a solution on X = Y 2 

line and one solution for X > Y 2 , so n + 2n = 3n cusps on lens plane. 
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(iii) \QP\<\Qn\<\Q®\ 

After butterfly singularity appears, there is a solution on X = Y 2 line 

and two solutions for X > Y 2 , so n + 2 x In = 5n cusps on lens plane. 

(iv) \Q®\<\Qn\ 

After swallowtail singularity appears, there is a solution on X = Y 2 
line and no solution for X > Y 2 , so n cusps on lens plane. 

(Ill) for n > 4 

« IQn|<|Q( 0) | 

There are three solutions on X = Y 2 line and one solution for X > Y 2 , 
so 3n + 2n = 5n cusps on lens plane. 

(ii) |Q(°)|<|Q n |<|Q«| 

After beak-to-beak singularity appears, there is a solution on X = Y 2 
line and one solution for X > Y 2 , so n + 2n = 3n cusps on lens plane. 

(iii) |Q«| < \Q n \ 

After butterfly singularity appears, there is a solution on X = Y 2 line 
and no solution for X > Y 2 , so n cusps on lens plane. 

5 Two Point-Mass Model 

Critical curve produced by the lens galaxy of PG1115+080 makes us imagine 
that the lens galaxy may be approximated by the two point-mass lens [Fig. 3]. 



Fig.3 



This model is interesting as its own right from the following reason. In the 
previous section we studied the complex lens object by decomposing it into 2 n - 
pole contribution. Two point-mass model is very simple one but is composed of 
composite multipoles. 

So in this section we analyse this model in the analogous manner as in the 
previous section. It is found that this model is qualitatively different from the 
MPE model. 

This model has been discussed by Erdl and Schneider [[J. A straightforward 
application of catastrophe theory to this model still enforces numerical calcula- 
tions. However we can derive singularities in simple analytical ways by pursuing 
the critical conditions of f(r,p) defined in Eq. (|5.8|) . 

The surface mass density of two point- mass lens is given by 

£(£) = M[fi A 5 2 (x-a) +p B 5 2 (x + a)] , (5.1) 



here M is total mass and \xa + A*b — 1- o and —a are locations of two point 
masses in the lens surface. 

Fermat potential </> is given by 

~ ! 777- 

(p = ^\ z - z s \ 2 - ^(ha\v-\z - x\ 2 + Hb\v-\z + x\ 2 ) , (5-2) 

in complex representation with m = -|r- and x — Q-i + w&2- 
We can rescale the potential into 

<j>= — = -\u-v\ 2 - -[p A ln|w - £| 2 + fi B \n\u + £| 2 ] (5.3) 

up to irrelevant constant. Here u = -£=, t> = -r= and £ = -£=. 
Lens equation ^ = gives 

U = U-( 7 + —T • 5 - 4 

Critical curve which is given by D = | — |0 UU | 2 = in general system of 
coordinates takes the following form, 



(u-0 2 (« + 2 



Introducing the new complex variable U = ^p-, £ in Eq.( p.5| ) is factorized out 



as 



2 |t/ 2 -4// B t/ + 4// B | 
lel = |^-2)P • (5 - 6) 



Using the polar coordinates U = re , critical curve Eq. (|5.6|) is rewritten as 

|£| 4 = /(r,cosfl), (5.7) 

where 



_ 16/i B r 2 p 2 - 8p B r(r 2 + 4// B )p + r 4 + 8/i B (2/i B - l)r 2 + 16/4 
/lr ' PJ " r 4 (r 2 -4rp + 4) 2 l 



with p = cos#. 

londition |g = 1 

5/ 



Beak-to-beak condition |^ = becomes 



«, (5.9) 



from Eq.( |5.7| ). Eg. (|5~^ ) is equivalent to 



Re(— )= d ^ I X - p2 — - (5 10) 

<9f/ <9r r <9p 
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and 



and, therefore, to 



for sin 9 7^ or to 



, ,df . df pdf 

Im(-) = sm£(----) = 0. (5.11) 



case 1) — - = tt - = (5.12) 

or op 



df 
case ii) sin9 = — — = . (5.13) 

or 

Beak-to-beak appears at the points where four real parameters r, p, p, B and 
|£| satisfy three equations (|5.71) and (|5.12|) or (|5.13|) . 



Thus metamorphoses of the caustics are appeared in the critical behaviours 
of f(r,p). 

Firstly we consider case i). 

df_ _ = 8(1 - p B )(16p B p - 12p B r + r 3 ) 
dp r 2 (r 2 — 4rp + 4) 3 

From the definitions, 1 — p B > 0, r 2 — 4rp + 4 > (|p| < 1 for case i)) and / 
takes the minimum value F m (r) at 

3 r 3 

p = -r = m(r) (5.15) 

4 16pB 



for fixed r. 



p B (r 2 -Ap B ) 2 



F m (r) ee /(r, m(r)) = ^-f^- 2 7 (> 0) . (5.16) 



We proceed to study the critical behaviour of F m (r). 

dF m (r) 4p B (r 2 -4p B ) 



Or r 5^ r 4 _ 8^^7-2 _|_ igp 2 ^) 

where 



g(r 2 ) , (5.17) 



ff ( r 2) = r e _ i2p B r 4 + 48p|r 2 - 64p| . (5.18) 

g(r 2 ) is a monotonically increasing function of r 2 and the root of g(r 2 ) = is 
given by r 2 , 

r? = 4/4 /3 (pf + pf) = 4p 2 /5 . (5.19) 

The behaviour of F m (r) is depicted in Fig. 4. 
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Fig.4 



Nextly we consider case ii) where p = ±1. We define f± = f(r, ±1). 
Then 

= (r 2 + 4/i B r + 4/xb) 2 
■>- ~ r 4( r + 2)4 

is a monotonically decreasing function of r and has no critical point. On the 
other hand, 

_ (r 2 - Ap B r + 4/i^) 2 

/+ " r 4 (r _ 2)4 ^.20) 



and 

M, A 

r 2 - Ap B r + Ap B ) /i(r) , (5.21) 



df+ 4 , 2 



dr r 5 (r — 2) 5 

where 

h(r) = r 3 - 6/i B r 2 + Ylp B r - 8p B . (5.22) 

h(r) is a monotonically increasing function of r and has only one real solution 
r 2 to h(r) = 0, 



r 2 = 2p L ^(p^ + /4 /d - vTv-b*) = —f~ • (5-23) 



1/3/ 2/3 2/3 1/3 1/3 \ _ 2/X s 



The behaviour of /+ is depicted in Fig. 5. 



Fig.5 

The above discussions in case i) and case ii) together with the facts that 

f(0,p) = oo , /(oo,p) = and (5.24) 

/(2,1) = oo 

are summarized as the contour map of f(r,p) [Fig. 6]. 

Fig. 6 
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It is found that the critical curve changes its topology when we cross the 
contour L\ (L 2 ) passing through C\ (C 2 ). We denote the £ values corresponding 
to the contour L\ (L 2 ) as £1 (£2). Then 



and 



N 4 = f\ch = 7^1 (5-25) 



I6| 4 = /Ic 2 = ^, (5.26) 



which agree with the results of Erdl and Schneider || . 

6 The Application to the Multiple Quasar PG1115+080 

In this section, we apply our MPE model to the lensing galaxy of multiple quasar 
PG1115+080. We use deflection potential up to the 2 3 -pole expansion here. 

About this lensed quasar, we have got the observed data of positions and 
relative amplifications of the images 0. (We have adopted the new data for only 
Z\ from Angonim-Willaime et al. [ID]) The positions of the images are as follows; 



0(A 1 ) = (-0".94, -0".73) 

9{A 2 ) = (-l".ll,-0".27) 

9(B) = (0".72, -0".60) (6.1) 

6(C) = (0".33, 1".35) 

And the relative amplifications of the images are as follows; 

IM^iVMC)! = 3.22 

IM40MC9I = 2.49 (6.2) 

HB)/f,(C)\ = 0.64 

Lens equation up to the 2 3 -pole expansion is given by 

1 D 1Q 3T . . 

v* -u* + - + — + -1 + — = 0. 6.3 

u u z u a vr 

Here v is source position, u is image position, and D, Q, T are dipole, 
quadrupole and 2 3 -pole, respectively. These are all dimensionless complex num- 
bered quantities normalized by suitable power of mass (See Eq.( |3.7| )). Image 
amplification is given by 
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.„ ,m 2d 6g 12i l2l , ,„ .. 

c = n-i F + F + ? + ii-rr- (6.4) 

We obtain values of the parameters included in Eg . ( |Q|) which reproduce the 
observed data well. These parameter values are as follows; 



source position = (0". 02785, 0".01507) 
m = 3.496 x KT 11 [rad 2 ] 
\D\ = 0.1824 , 9 D = 88°.81 (6.5) 

\Q\ = 0.08724 , 9 Q = 130°. 6 

h 



\T\ = 3.407 x 10~ 3 , 9 T = 115°. 5 



If we assume Einstein-de Sitter universe, this mass would be M = nT> cr Dfm = 
1.483 ± 0.024 h~ x x 10 11 M . Here redshifts of the lensing galaxy and the source 
quasar are z\ = 0.294 ± 0.005 and z s = 1.722, respectively. Hubble constant is 
H = 100 h km/sec/Mpc. 

Calculated image positions are 



9{A 1 ) = (-0".9383,-0".7317) 

9(A 2 ) = (-1". 110,-0". 2748) 

0(B) = (0".7205, -0".6004) (6.6) 

9(C) = (0".3294, 1".350) 

9(D) = (0".5915, -0".4450) 

9(E) = (-0".4494, // .08641) . 



And calculated relative amplifications are 



\n{A 1 )/n(C)\ - 


= 3.113 


\V(A 2 )I^(C)\ - 


= 2.892 


HB)/H(C)\ - 


= 0.7467 


\n{D)ln(C)\ - 


= 0.1780 


HE)/KC)\ - 


= 0.005871 . 



(6.7) 



Then, we should check the validity of approximation in multipole expansion 
up to 2 3 -pole. It must hold the following conditions. 

i » \£l » M » m ( 6.8) 

\u\ \ur \u\ 6 



And the best fit parameters ( |6.6| ) give following values. 
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\D\ \Q\ \T\ 9 . . 

l —t » 0.2609 , f^ « 0.2103 , f-1 « 1.621 x 10~ 2 6.9 

|it| \u\ 2 \u\ 6 

This result roughly holds the condition P& ^$> -hX. So this MPE model up to 
2 3 -pole has been proved to be valid up to 2 3 -pole moment. 

We performed parameter fitting of this lensed quasar using MPE model up to 
quadrupole in the previous article J7| . The values of source position, dipole and 
quadrupole in Eq. (|6.6|) are similar to those values in previous fitting. So they are 
consistent. 



7 The Relations with Other Models 

We have shown in the previous sections that deflection potential up to the 
quadrupole expansion reproduces multiple images well. Its critical curve takes 
the similar form produced by two point-mass lens ||. 

So it may be meaningful to check that two point-mass lens model on the single 
lens plane well resembles with the lens of PG1115+080. 

Also it is very helpful for the study of lensing object that we analyze the 
same lensing system by qualitatively different models. As a typical example, we 
consider softened power-law spheroid models. 

7.1 Two point-mass lens 

Fermat potential given by Eq. (|5.2|) is expanded by multipoles as 



Here 



I — * — *| 9 I — ♦) LLj^Ju (J'ijJb Jb 

— \x — y\ — mmx — -— — 

2 \xr \x 



x — y\ -mm|s|-|^- —^ ... . (7.1) 



M 



m = — — — 

7TZj cr 




di = m^2 ^kdki 

k 


{k = A,B) 


lij - m Yl ^k{d k id k j ■ 

k 


141% , 

- 2 Sy) 



(7.2) 



and so on, which are the real representations of Eqs.( |3.3j ) 

In this case we have six independent parameters M , \ia , ^a ,ds besides source 

positions y. Here we assume that MPE model up to quadrupole moments is good 

approximation as we have shown and that two point-mass lens can be considered 

to be their prototype. Then we must be able to fit the parameters M , jUa , ^a ,d>B 
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from Eq.( |7.2| ) so as to give the same m ,di , qij as the MPE model. This is always 
possible. For Eq. (|7.2j ) consists of five equations and we have six parameters to 
be adjusted. 

Whether two point-mass model is good approximation or not, therefore, de- 
pends on the fact that the higher multipole terms subsequent to quadrupole in 
Eq. (|7.1|) is small or not. Here higher multipole terms is fixed by M , fj,^ , d& , ds 
determined by Eq.( |7.2|) . 

As is easily checked it is impossible. Therefore two point-mass lens model 
cannot be a good approximation of PG1115+080. 

7.2 Softened power-law spheroids (SPLS) 

This model is also called elliptical lens model whose mass density is given by 

p(a)=p (l + ^)- 1 (7.3) 



c 



with 



x' 2 + y' 2 + ~^— 2 = a 2 (7.4) 

1 — e z 

The principal axes of spheroid (x 1 , y',z') are related with the coordinate of 
optical axis and lens surface (x,y,z) by 

/ 1 \ / x \ 

cos 7 —sin 7 j y . (7.5) 

\ sin 7 cos 7 ) \ z ) 

That is, optical axis z is tilted to z' by 7. 

Here r c is the core radius and po is the constant central density, e is the 
eccentricity. For k < 3 total mass is divergent, mass density has a cut-off 

p(a) =0 for — > n . (7.6) 




This model has been applied to PG1115+080 by Narasimha et al. ]TTJ. How 



ever they used the observational data different from ours, therefore we cannot 
simply compare their results with ours. Elliptical lens model is qualitatively 
different model from MPE model. The former is transparent and the latter non- 
transparent. If we are only satisfied by reproducing image positions, both model 
may pass our requirements. So we must try parameter fitting to more various 
and precise data. 

The detail of this model has been discussed in a separate form 0, so we just 
quote the result [Table 1, Table 2]. 
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Table 1 



Table 2 



SPLS has a large fitting zone for its transparency compared with MPE. 

As a whole this model gives rather non-compact lensing object unlike the 



compact aspect of Narasimha et al. [O. This seems to be incompatible with the 



observed lensing object fT2| . Observed image amplifications are better fitted by 



MPE model than by elliptical lens model. 

8 Discussions 

MPE model is not so popular as the elliptical lens model with or without external 
shear since most of the gravitational lensings are thought to be caused by the 
transparent global systems composed of many galaxies. However, the recent 
observations support that the lens of PG1115+080 is a compact object whose 
optical extent is within the observed innermost image B ~ lOh^ 1 kpc |J [ [12] . 
MPE model is neither transparent nor nonsingular. Therefore it is free from the 
general theorems about gravitational lensing (See chapter 5.4 of Schneider et al. 
[|TJ). Indeed, in Fig. 3 we have 4 images of odd parity and 2 images of even parity. 
There gives rise, of course, no inconsistency. 

We must wait further refinement of model check in order that our model 
parameters are physically affirmative. For that purpose, we need more precise 



observations such as image shape, time delays and spectroscopy profile [13|], etc. 
In this connection it is in order to comment on time delay. 
Time delay is given by 



t= -(1 + ^)^0 • (8-1) 

c Di s 



For Qn = 1 case 



2C (1 + Zj )y/1 + Zi - (1 + Zi)^/T+ 



Zj 



A, = ir K 3 '\: * , 7 V J ■ (8-2) 



For fin = case, 



H (l + z i ) 2 (l + * 

Zi < Zj , i,j = source, lens, observer ) 
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Time delays for each case are given in Table 3. 



Table 3 



For comparison we have also listed the respective time delay based on the 
elliptical lens model [Table 2]. Estimations by MPE model are roughly coinci- 
dent with the observation by Vanderriest et al.[Q- On the contrary, estimations 
by SPLS in Table 2 are too short, which denies the possibility that its amplifi- 
cation discrepancy with observation might be remedied by the time variation of 
amplification of source QSO. Naively we may conclude that the lensing object of 
PG1115+080 is compact also from the theoretical side. 

Consequently MPE model deserves serious surveys. In this connection it 
should be remarked that the dipole term plays an important role in our theory : 
The presence of dipole moment means that the origin of the coordinates (which 
is the optical center) is away from the center of mass of the lensing galaxy. 
There are the arguments which insist that this deviation may be the cause of 
spiral arms. Also MPE model gives better fit to the observational data than 
elliptical lens model in many points ||. This may supports the view point that 
the lensing galaxy is a spiral galaxy, whose equipotential contour is different from 
spheroid. There is the observational support of this viewpoint from the analysis 



of spectroscopy |I3|p5|. From these facts, non zero dipole component may play 
any role in the formation of spiral arms. 

Though we applied it only to PG1115+080 in this article, MPE model is 
very useful for the wide class of lensing phenomena. It gives systematically the 
multipole moments of the lensing object, irrespective to isolated lens or not. This 
serves crucially to the construction of more concrete and more precise structure 
of lensing object. 
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Figure Captions 

Fig.l The region that satisfies the cusp condition. Solid line is X = Y 2 , dotted 
line is Re(cusp)=0 and dashed line is lm(cusp)=0. Three lines pass through 
the points P\ and P2. The cusp condition is satisfied in the region where 
the dashed line is in the common region of the interiors of the solid and the 
dotted line. Between Pi and P2, the dashed line passes very close to but in 
the exterior of the dotted line. 

Fig. 2 The behaviours of the critical condition. Solid line is X = Y 2 , dashed 
line is lm(cusp)=0 and dotted line is the critical condition. Fig.2(a) is the 
n = 2 case. The dotted lines of (0), (1) and (2) correspond to the cases (0), 
(1) and (2) discussed in the text, respectively. Fig. 2(b) is the n = 4 case. 

Fig. 3 Shapes of images, critical line and caustics for the multiple quasar PG1115+080 
(Kakigi et al. 1995). Crosses (dotted regions) are the observed (calculated) 
image positions. It has been assumed that source QSO is spherical. 

Fig. 4 The critical behaviour of F m (r). Critical points Cq and C\ in (r, p) plane 
are (r ,p ) = (2/i^ /2 ,/i^ /2 ) and (r^p^ = (2/4/ 3 <5 1/2 , ^(2p^ 3 - fi 1 / 3 )), 
respectively. We can set ps > Ha without spoiling generality and therefore 
1 > pi > 0. 

Fig. 5 The critical behaviour of /+. Critical point C2 is given by C2 = (r 2 , P2) = 

Fig. 6 The contour map of f(r,p). Origin is the point of heavier point mass. 
f(r,p) diminishes in the direction of =K C\ and C2 are the points where 
beak-to-beaks appear. 

Table 1 Fitting parameters of SPLS ; (a) k=3 case and (b) k=2 case. 

Table 2 Calculated image positions (and time delays lagging behind C in h~ l 
days) of SPLS ; (a) k=3 case and (b) k=2 case. 

Table 3 Time delays lagging behind C in h" 1 days of MPE model. 
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esiii7 


n 


r c 


/toVl — e 2 


Total Mass 


Source Position 


0.4 


5 


6.13 h^ 1 kpc 


2.58xl0~ 23 h 2 g/cm 3 


1.48xlO i2 M 


(0".024, -0".037) 


0.4 


10 


6.07 h- 1 kpc 


2.56xl0~ 23 h 2 g/cm 3 


2.13xlO 12 M 


(0".024, -0".037) 


0.5 


5 


4.78 h- 1 kpc 


3.48xl0~ 23 h 2 g/cm 3 


9.43xlO n M 


(0".038, -0".058) 


0.5 


10 


4.71 hr 1 kpc 


3.50xl0- 23 h 2 g/cm 3 


1.37xlO 12 M 


(0".039, -0".058) 


0.6 


5 


3.84 h- 1 kpc 


4.66xl0~ 23 h 2 g/cm 3 


6.56xlO n M 


(0".057, -0".085) 


0.6 


10 


3.75 h~ x kpc 


4.74xl0~ 23 h 2 g/cm 3 


9.3OxlO n M 


(0".057, -0".084) 
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n 


r c 


KoVl — e 2 


Total Mass 


Source Position 


0.3 


5 


6.00 h- 1 kpc 


1.79xl0~ 23 h 2 g/cm 3 


2.62xlO 12 M 


(0".014, -0".020) 


0.3 


10 


5.46 h~ l kpc 


1.87xl0" 23 h 2 g/cm 3 


4.83xlO 12 M 


(0".014, -0".021) 


0.4 


5 


4.23 h- 1 kpc 


2.71 xl0~ 23 h 2 g/cm 3 


1.39xlO 12 M 


(0".025, -0".036) 


0.4 


10 


3.97 h~ l kpc 


2.71 xl0~ 23 fi 2 g/cm 3 


2.7OxlO i2 M 


(0".025, - 0".035) 


0.5 


5 


3.13 h- 1 kpc 


3.93xl0~ 23 h 2 g/cm 3 


8.15xlO n M 


(0".039, -0".055) 


0.5 


10 


2.84 h- 1 kpc 


4.13xl0~ 23 h 2 g/cm 3 


1.51xlO 12 M 


(0".040, -0".054) 
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n 


A x 


A 2 


S 


C 


0.4 


5 


(0".51, 1".08) (4.0) 


(0".95, 0".63) (4.1) 


(-0".88, 0".31) (6.0) 


(0".23, -1".37) 


0.4 


10 


(0".51, 1".08) (3.9) 


(0".96, 0".62) (4.0) 


(-0".89, 0".31) (5.9) 


(0".24, - 1".37) 


0.5 


5 


(0".53, 1".06) (6.1) 


(0".97, 0".61) (6.3) 


(-0".89, 0".30) (9.4) 


(0".24, - 1".37) 


0.5 


10 


(0".53, 1".06) (6.2) 


(0".97, 0".61) (6.3) 


(-0".89, 0".30) (9.6) 


(0".25, - 1".37) 


0.6 


5 


(0".53, 1".05) (9.0) 


(1".00, 0".57) (9.2) 


(-0".89, 0".28) (13.9) 


(0".26, - 1".37) 


0.6 


10 


(0".53, 1".05) (8.8) 


(1".00, 0".57) (9.1) 


(-0".89, 0".28) (13.8) 


(0".26, - 1".37) 


Table 2 (a) 
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n 


A, 


^2 


B 


C 


0.3 


5 


(0".52, 1".08) (2.1) 


(0".94, 0".64) (2.1) 


(-0".88, 0".32) (3.2) 


(0".23, -1".37) 


0.3 


10 


(0".56, 1".04) (2.2) 


(0".92, 0".67) (2.3) 


(-0".90, 0".31) (3.4) 


(0".23, -1".35) 


0.4 


5 


(0".54, 1".06) (3.9) 


(0".96, 0".62) (4.0) 


(-0".88, 0".31) (6.1) 


(0".26, - 1".36) 


0.4 


10 


(0".52, 1".07) (3.8) 


(0".92, 0".61) (4.0) 


(-0".89, 0".30) (6.1) 


(0".26, -1".36) 


0.5 


5 


(0".54, 1".05) (5.8) 


(0".99, 0".58) (6.0) 


(-0".89, 0".29) (9.3) 


(0".28, -1".36) 


0.5 


10 


(0".54, 1".05) (5.8) 


(0".99, 0".57) (6.0) 


(-0".89, 0".29) (9.4) 


(0".28, -1".36) 



Table 2 (b) 



A (h- 1 days) 


C-A 1 


C-A 2 


C-B 


C-D 


C-E 


Q = l 


9.757 ±0.200 


11.448 ±0.234 


19.407 ±0.396 


19.225 ±0.392 


3.076 ±0.063 


tt = 


10.540™ 


1Z.OO/_ 264 


20.964+^47 


20.767™ 


o qoq+U.UV2 
°- ozo -0.071 
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